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With growing datasets collected by the CDF II experiment, studies of the spectroscopy of mesons 
containing heavy quarks become more exciting. The CDF experiment has good capabilities in 
both charm and bottom sector. This capability allowed also to contribute to the study of the 
zoo of states called X,Y,Z. In this area we present a recent update of the mass measurement of 
X(3872). The result m(X(3872)) = 3871.61 ±0.16 ±0.19MeV/c 2 is currently the most precise 
measurement in the world. In addition, we report evidence for a new narrow resonance, 5 / (4140), 
the first to be seen in the J/\j/^> decay mode, using 2.7fb _1 of exclusive B + — > J/\j/(j)K + decays. 



European Physical Society Europhysics Conference on High Energy Physics 
July 16-22, 2009 
Krakow, Poland 



* Speaker. 

"'"On behalf of the CDF Collaboration. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos . sissa . it / 



Heavy quark meson spectroscopy at CDF 



Felix WICK 



1. Introduction 

In recent years, several states with charmonium-like decays were discovered which do not 
fit properly into the established charmonium picture. These states, called X,Y,Z, are candidates 
for exotic mesons beyond the conventional quark-antiquark-model (qq). Possible interpretations 
are quark-gluon-hybrids (qqg), four-quark states (qqqq), molecular states composed of two usual 
mesons, or glueballs. In the following, two CDF analyses in this context are considered, a measure- 
ment of the X(3872) mass in J/yn + n~ decays [JTJ] and an evidence for a narrow near-threshold 
structure in the J /yip mass spectrum, called 7(4140), using B + — > J / \j/(j)K + decays [0]. 

Although the existence of X(3872) is well established, its nature is still unclear. The two 
most popular models are a molecular state composed of D° and D*° mesons and a four-quark state 
(ccqq). The molecular model is motivated by the closeness of the measured mass to the sum of 
the D° and D*° masses and could be rejected if m(X(3872)) was found larger than this sum. The 
four-quark state hypothesis predicts the existence of two distinct particles with different light-quark 
content and slightly different masses. Because of these predictions, an experimental examination 
of the X(3872) mass and its structure allows to test the two hypotheses. 

Consisting of two vector mesons, the final state J/y<j> is a good channel to search for an 
exotic meson. Because the invariant mass in this channel is high enough for open charm decays, 
a charmonium state is very unlikely to be observed. Using exclusive B + decays to J/ytpK + , a 
strong background reduction can be achieved by exploiting the long B meson lifetime and the 
identification of the additional kaon. 

The data used for the reported analyses were collected by the CDF II experiment at the Teva- 
tron proton-antiproton collider at Fermilab. CDF II is a multipurpose detector with an excellent 
tracking and mass resolution. Therefore it is well suited to perform the described measurements. 

2. X(3872) Mass 

X(3872) candidates are reconstructed in the decay channel J /yii + K~ with J/y — > }JL + pL~ . 
y(2S) decays to the same final state and is used as a control sample. In order to separate signal from 
background, a neural network is trained by means of simulated X(3872) events and the sidebands 
in the data mass spectrum. The most important input quantities for the neural network are the Q 
value of the decay, the transverse pion momenta, the kinematic fit quality of the X (3872) candidate 
and muon identification quantities. Subsequently, cuts on the output of this neural network and the 
number of candidates per event are performed which select around 34500 y(2S) and 6000 X (3872) 
signal events. 

In order to test whether there is evidence for more than one state, a binned maximum like- 
lihood fit to the mass spectrum is performed, where the signal is described by a nonrelativistic 
Breit-Wigner function convolved with a Gaussian resolution and the combinatorial background is 
modeled by a second order polynomial. Thereby the resolution is obtained from MC simulations 
and the Breit-Wigner width is fixed to 1.34MeV/c 2 (average of [Q| and [Q]). But additionally, the 
overall width of the signal shape is scaled by a factor t which serves as a free parameter in the 
fit. The parameter t is also determined for ensembles of simulated experiments consisting of two 
X(3872) states with given mass differences which are generated using the same fit model as in data 
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Figure 1: (a) Distribution of the width scale t from pseudoexperiments assuming a single X(3872) state. 
Also shown is the value measured in data, (b) Upper limits on the mass difference between two possible 
states in dependence of the low-mass signal fraction. 
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X(3872) Mass Measurements in J/\|f7nr. Mode 



Belle 

3871 .46 ± 0.37 ± 0.07 MeV/c 2 
BaBar (B + ) 

3871 .30 ± 0.60 ±0.10 MeV/c 2 

BaBar (B°) I 

3868.60 ± 1 .20 ± 0.20 MeV/c 2 
DO 



3871.80±3.10±3.00 MeV/c^ 
CDF Old 

3871 .30 ± 0.70 ± 0.40 MeV/c 2 
CDF new (preliminary) 
3871.61±0.16±0.19MeV7c 2 



average w/o new CDF result 
3871. 24 ± 0.29 MeV/c ? 
average with new CDF result 
3871 .46 ± 0.19 MeV/c 2 



m(D D )+m(D*°) [PDG] 
3871.81 ±0.36 MevVc* 



b) 



3868 3869 3870 3871 
X(3872) Mass ( MeV/c 2 ) 



Figure 2: (a) Invariant mass distribution of the X(3872) candidates with the fitted function (solid red line) 
and the background model (dashed line), (b) X(3872) mass measurements in the decay to J / \j/7t + n~ . 



but varying several quantities, like the Breit-Wigner width, according to their uncertainties. Figure 
[I] (a) shows the comparison of the measured value t with its distribution resulting from an ensemble 
of pseudoexperiments assuming a single state (Am = OMeV/c 2 ). Here no evidence for two states 
can be found and an upper limit on their mass difference can be set by finding the value for which 
the scaling parameter t of 90% (95%) of the pseudoexperiments is higher than the measured t. In 
figure [I] (b) this limit is shown in dependence of the low-mass signal fraction f\. For f\ = 0.5 the 
limits are m < 3.2MeV/c 2 and m < 3.6MeV/c 2 at 90% and 95% confidence levels, respectively. 
This can be compared to the prediction of Am = (8 ± 3) MeV/c 2 by the four-quark state model of 
Maiani et al.[Q]. 

Since the data are fully consistent with a single state, the mass of the X (3872) can be measured 
by performing an unbinned maximum likelihood fit using the same fit model as described above. 
The result of this fit is shown in figure || (a). The momentum scale uncertainty as main source of 
systematic errors can be estimated by means of the difference between the measured y(2S) mass 
and the corresponding world average value. As can be seen in figure |2| (b), the resulting value 
to(X(3872)) = 3871.61 ±0.16 ±0.19MeV/c 2 is the most precise single measurement to date and 
lies below m(D°) J r m(D*°). So the molecular model is still possible. 
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Figure 3: Invariant mass distributions of (a) J/\j/<pK + and (b) K + K with the fitted functions (solid blue 
lines) as described in the text. 



3. Evidence for Y (4 140) 

In order to build B + -^J/\j/(j)K + candidates, first J/y — > and — > K + K candidates are 

reconstructed which are then combined with an additional charged track with kaon mass hypothe- 
sis. Thereby the reconstructed J j\\f and masses are required to lie within 50MeV/c 2 (J/y) and 
7MeV/c 2 (0) of the corresponding world average values. The combinatorial background can be 
reduced significantly by cutting on the B + decay length in the transverse plane (L xy (B + )) to exploit 
the long B meson lifetime. In addition, a kaon identification quantity can be used to obtain a further 
background reduction. For that purpose, the information about the energy loss dEldx in the drift 
chamber and the information from the Time-Of-Flight detector are summarized in a log-likelihood 
ratio LL/?Kaon- The cuts L XJ (B + ) > 500/^m and LLR^ aon > 0.2 are chosen by optimizing the quan- 
tity S/yJS + B, where S and B are the numbers of B + signal and background events, respectively. 
Figure ^] (a) shows the resulting J/\j/(f)K + invariant mass spectrum. A fit with a Gaussian signal 
and a linear background function yields 75 ± 10 signal events which is the largest sample to date. 
For the examination of the J/yftf) spectrum, only candidates within ±3a(17.7MeV/c 2 ) around the 
nominal B + mass are selected. Figure |3| (b) shows the B + sideband-subtracted — > K + K invariant 
mass spectrum without <p mass window requirement. The fit function is a P-wave relativistic Breit- 
Wigner convolved with a Gaussian to account for the detector resolution. As there is no significant 
background contribution, the B + — > J/\j/K + K~K + final state is well described as J/\j/<pK + . 

In a Dalitz Plot of the considered decay (see figure |] (a)), an enhancement near the thresh- 
old of m 2 (J/\j/^)) can be seen. Figure || (b) shows the spectrum of the mass difference AM = 
m(pL + ii~K + K~) — m{pL + pL~) using 73 events with AM < 1.56GeV/c 2 . An unbinned maximum 
likelihood fit is performed, where the enhancement is described by the convolution of an S-wave 
relativistic Breit-Wigner function with a Gaussian resolution with the rms fixed to 1.7MeV/c 2 ob- 
tained from MC simulations, and the background is modeled by three-body phase space. The fit 
yields 14 ±5 signal events and returns a mass difference AM = (1046.3 ±2.9) MeV/c 2 and a width 
of (11.7+^)MeV/c 2 . 

In order to estimate the probability of a creation of such a signal due to background fluctu- 
ations, simulations of the background distribution are performed and the number of trials which 
produce a signal with a log-likelihood ratio — 21n( = Sfo/-Sfmav) of the null hypothesis fit and the sig- 
nal hypothesis fit larger than the measured value are counted. Thereby the mass can be anywhere in 
the considered AM window and the width has to be larger than the detector resolution and smaller 
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a ) m 2 (J/y<W (GeV 2 /c 4 ) b) AM (GeV/c 2 ) 

Figure 4: (a) Dalitz plot of the final state J/\j/^)K + in the B + mass window. The boundary shows the 
kinematically allowed region, (b) Distribution of the invariant mass difference AM = m{pL + pL~ K + K~) — 
m{ii + pL~) in the B + mass window with the fitted function (solid red line) and the background model (dashed 
blue line). 

than ten times the observed width. In addition, the log-likelihood ratio in data is decreased by 
modeling the combinatorial background in the B + mass window separately as flat spectrum in AM 
instead of using pure three-body phase space background. This procedure leads to a significance 
of the enhancement of about 3.8a. 

With the world-average value for the J/y mass, the mass and width of the observed structure 
are (4143.0 ± 2.9 ± 1.2) MeV/c 2 and (11.7+|^±3.7)MeV/c 2 , respectively, where the systematic 
uncertainties are estimated by varying the fit model. 

4. Summary 

The presented measurement of the X(3872) mass in its decay to J/\j/7t + 7i~ with a value of 
(3871.61 ±0.16(stat) ±0. 19 (syst)) MeV/c 2 is the most precise measurement to date. Additionally, 
an upper limit on the mass difference of two hypothetical X(3872) states is derived. 

In the J I mass spectrum of exclusively reconstructed B + — > J/y<pK + decays, evidence for 
a narrow near-threshold structure with a significance in excess of 3.8a is found. The measure- 
ments of the mass and width of this structure yield (4143.0 ±2.9(stat) ± 1.2(jyjf))MeV/c and 
(U.7± % 5 l(stat) ± 3.7(yytf))MeV/c 2 , respectively. 
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